Linear quantum trajectories: Applications to continuous projection measurements 
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We present a method for obtaining evolution operators for linear quantum trajectories. We apply 
this to a number of physical examples of varying mathematical complexity, in which the quantum 
trajectories describe the continuous projection measurement of physical observables. Using this 
method we calculate the average conditional uncertainty for the measured observables, being a 
central quantity of interest in these measurement processes. 
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I. INTRODUCTION 

Quantum master equations, which govern the evolu- 
tion of a density matrix representing the state of a phys- 
ical system, have a wide application in quantum dissipa- 
tion and continuous measurement theory [H, [2, 0] . They 
describe the evolution of a quantum system that is in- 
teracting with an environment which, due to the interac- 
tion, may absorb energy from the system (dissipation), 
and will continually provide information about the state 
of the system (continuous measurement). A classic exam- 
ple of a system interacting with an environment is that 
of a single mode of an optical cavity which is allowed 
to leak out of the cavity via an imperfect end-mirror. 
The photons in the cavity leak out over time, and these 
may be detected with a photo-detector. The environ- 
ment consists of the continuum of optical modes out- 
side the cavity, and provides a mechanism for dissipation 
and continuous measurement. A master equation would 
describe the evolution of the system averaged over all 
the possible times at which the photons may be detected 
leaving the cavity. However, the master equation may be 
rewritten in an equivalent form as a stochastic equation 
which describes the evolution of the system for each set 
of photo-detection times [i| . Each possible realisation of 
the stochastic equation corresponds to a set of detection 
times, or more generally, to a particular set of measure- 
ment results. Each set of results is termed a quantum 
trajectory and the stochastic equation is said to un- 
ravel the master equation. The kind of stochastic process 
appearing in the equation will depend upon the kind of 
measurement process. For photo-detection of the output 
of a cavity mode the stochastic process is a point pro- 
cess 16 1, while for homodyne detection it is a Wiener pro- 
cess [7|. However, the master equation, giving the overall 
average evolution, docs not depend upon the choice of 
measurement. In other words, there are many different 
ways to unravel any particular master equation. Here we 
will be concerned with stochastic equations which con- 
tain the Wiener process. 

The fact that quantum master equations may be 
rewritten as Linear Stochastic Equations (LSE's) for the 
quantum state vector, referred to alternatively as linear 
quantum trajectories, has been known in the mathemati- 
cal physics literature for some time Q , but has only fairly 



recently seen exposure in the physics literature 0, ITol . [TlT | , 
where it has been common to use non-linear stochastic 
equations [Hi EH ■ The advantage of writing master equa- 
tions as LSE's, rather than the more familiar non-linear 
version, is that in certain cases it has been found that 
explicit evolution operators corresponding to these equa- 
tions may be obtained in a straightforward manner. How- 
ever, as far as we are aware, the only method that has 
been used to obtain evolution operators for these equa- 
tions to date is to choose an initial state which allows 
the stochastic equation for the state to be written as a 
stochastic equation for an eigenvalue, or which simplifies 
the action of the evolution operator [3, [l(| • In this paper 
we present a more general method for obtaining explicit 
evolution operators for these equations which makes no 
reference to the initial state. Naturally the resulting evo- 
lution operators contain classical random variables. The 
complexity of the stochastic equations which govern these 
classical random variables depends upon the complexity 
of the commutation relations between the operators ap- 
pearing in the LSE. If the complexity of the commutation 
relations is sufficiently high then the stochastic equations 
governing the classical random variables become too com- 
plex to solve analytically. Nevertheless, even if this is 
the case, the form of the evolution operator provides in- 
formation regarding the type of states produced by the 
LSE, and the problem is reduced to integrating the clas- 
sical stochastic equations numerically. We also note that 
the solution to an LSE provides additional information 
to that contained in the solution to the equivalent mas- 
ter equation, because it gives the state of the system for 
each trajectory. For example, the variance of a system 
operator may be calculated for each final state (ie. for 
each trajectory), and this is referred to as the conditional 
variance as it is conditional upon the results of the mea- 
surement. The overall average of these variances may be 
then be calculated. The solution to the master equation 
allows us to calculate only the variance which is obtained 
by first averaging the final states over all trajectories, 
which is, in general, quite a different quantity. 

In the following we use as examples LSE's correspond- 
ing to the continuous measurement of physical observ- 
ables. A term of the form 



k[0, [0,p]} 



(1) 
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in a quantum master equation for the evolution of a den- 
sity matrix, p, for a quantum system S, describes a con- 
tinuous projection measurement of an observable O of S. 
The rate at which information is gained regarding the ob- 
servable is determined by k which is a positive constant. 
That a continuous measurement of a physical observable 
can be described in this way has been demonstrated by 
Barchielli and co-workers and also by Ueda et al. 
fl5| using a quite different approach. For the theory of 
continuous measurement the reader is referred to these 
works and references 0, [l?], El, 0- We refer to this 
measurement process as a continuous projection measure- 
ment because in the absence of any system evolution, the 
sole effect of this term is to reduce the off-diagonal ele- 
ments of the density matrix to zero in the eigenbasis of 
that observable. That is, it describes, in the long time 
limit, a projection onto one of the eigenstates of the ob- 
servable under observation. If, in addition, the observ- 
able commutes with the Hamiltonian describing the free 
evolution of the system under observation, then the free 
evolution does not interfere with this process of projec- 
tion, and the measurement is referred to as a continuous 
Quantum Non-Demolition (QND) measurement (3, [20j. 

Before we proceed we note the following points. The 
LSE which is equivalent to the general master equa- 
tion H2 

N 

P = - l j\H, P ] + J2(20 nP Oi - 0\fi nP - pOlO n ), (2) 

71=1 

where H is Hcrmitian and the O n are arbitrary operators, 
is 



variables. For an in-depth account of LSE's and their re- 
lationship to physical measurements we refer the reader 
to reference Q. 



II. OBTAINING EVOLUTION OPERATORS 
FOR LINEAR QUANTUM TRAJECTORIES 

A. General method 

We will explicitly treat here LSE's which contain 
only one stochastic increment. However it will be clear 
that this treatment may be easily extended for multiple 
stochastic increments. Let us write a general LSE with 
a single stochastic increment as 



d\ip) = (A dt + B dW{t))\ijj). 



(6) 



In this equation A and B are arbitrary operators. We will 
see that the complexity of the evolution operator will de- 
pend upon the complexity of the commutation relations 
between A and B. 

Let us define the integral of Wiener increments over a 
time At as AW(t). The probability density for AW(t) 
isi 



P{AW{t)) 



1 



V2ttAI 



-(AW(t)) 2 /(2At) 



(7) 



so that (AW(t)) = and ((AW(t)) 2 } = At. As a first 
step in obtaining an evolution operator for the LSE in 
Eq.([6]) we rewrite it in the form 



dm = 



N 



-xHdt - J2(OlO n dt - V20 n dW n {t)) 



\4>), 

(3) 

where the dW n (t) are independent stochastic Wiener in- 
crements which obey the Ito calculus relation dW n (t) 2 = 
dt 0. 

During evolution an initialy pure quantum state re- 
mains pure, but changes in a random way determined 
by the values taken by the Wiener process. The state 
at time t, \if>(t)) w , is not normalised, and the probability 
measure for the system to have evolved to that particular 
state at that time is given by [i[ 



(^W|VW)wdP w , 



(4) 



where dP w is the Wiener measure. That is, it is the joint 
probability measure for all the random variables that ap- 
pear in the expression for \ip(t)) w . It follows therefore 
that moments of system operators calculated with the 
equivalent master equation at time t are given by the 
expression 



(o) = / m)\o\m)^p^ 



(5) 



where O is the system operator in question, and dP w rep- 
resents integration over all possible values of the random 



\ip{t + dt)) 



e (A-(B 2 /2))dt e BdW(t), 
e Adt e BdW{t) U) 



(8) 



were we have defined A = A — B 2 /2. It is easily veri- 
fied that this is correct to first order by expanding the 
exponentials to second order and using the Ito calculus 
relation dW(t) 2 = dt. To first order the state at time 
t + At is therefore 



\^){t + At)) = e 



_ AAt BAW(t) 



\m), 



(9) 



so that the state at time t may be written 



lim 

At^O 



N 

IK 

' n=l 



e AAt & BAW„ 



where 



AW n 



At 

dW(t), 

(ra-l)At 



) W0)>, (10) 



(11) 



and N — > oo as At — > so that NAt = t is always true. 
To complete the derivation of the evolution operator we 
must take the limit in Eq. (|10[) . To do this we must com- 
bine the arguments of the exponentials which appear in 
the product, so that we may sum the infinitcssimals. We 
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will choose to do this by first repeatedly swapping the 
order of the exponentials containing the operator A with 
those containing the operator B. The simplest case oc- 
curs when A and B commute so that the problem essen- 
tially reduces to the single variable case, and we treat this 
in Sec lIIBI The simplest non-trivial case occurs when the 
commutator [A, B], while non-zero, commutes with both 
A and B, and we treat this in Scc lII CI In the final part 
of this section we examine a more complicated example 
in which the commutator [A, B] does not commute with 
either A or B. 



B. A QND measurement of photon number 

The mathematically trivial case occurs when A and B 
commute. A non-trivial physical example to which this 
corresponds is a QND measurement of the photon num- 
ber of a single cavity mode. Denoting the annihilation 
operator describing the mode by a, the free cavity field 
Hamiltonian is given by 0] 



H = hw{a)a + i), 



(12) 



in which u> is the frequency of the cavity mode, and the 
observable to be measured is O — a^a. With this we have 



A = -iw(a f a+i) - 2/c(a f a) 2 , 



B 



2kah 



(13) 
(14) 



in which k is the measurement constant introduced in 
Eq.([T]). As A and B commute the exponentials in Eq.lflO]) 
combine trivially and we obtain 



AN At , 



IV ) (*))w = lirn e 

At^O 



e At e BW(t) 



exp 
|V(0)> 



\m) 



(15) 



As the Wiener process, W(t), is a sum of independent 
Gaussian distributed random variables, W n , it is nat- 
urally Gaussian distributed, the mean and variance of 
W(t) being zero and t respectively. In a particular re- 
alisation of the stochastic equation Eq.©, the Wiener 
process will have a particular value at each time t, and 
as we mentioned above, the set of all these values corre- 
sponds to the trajectory that is taken by that particular 
realisation. The fact that to obtain the state at time t 
we require only the value of the Wiener process at that 
time means that we do not require all the trajectory in- 
formation, but just a single variable associated with that 
trajectory. For more complicated cases, in which the op- 
erators do not commute, we will find that other variables 
associated with the trajectory appear in the evolution 
operator. 

As the situation we consider here is a QND measure- 
ment, the phase uncertainty introduced by the measure- 
ment of photon number does not feed back to affect the 



measurement, so that the result is simply to decrease 
continuously the uncertainty in photon number, and the 
state of the system as t tends to infinity tends to a num- 
ber state. If we denote the evolution operator derived 
in Eq.(|15p by V(t), and start the system in an arbitrary 
mixed state p(0), then at time t the normalised state of 
the system may be written 



P(t) = 



V(t)p{Q)V\t) 
Tr{V(t)p(0)VHt)}' 



(16) 



As V(t) is diagonal in the photon number basis, we only 
require the diagonal elements of the initial density ma- 
trix to calculate moments of the photon number operator. 
Denoting the diagonal elements of the initial density ma- 
trix by p n , and the diagonal elements of V{t)V^{t) by V n , 
the variance of the photon number, for a given trajectory, 
is given by 



2 (f) = E„ n 2 p n V n ^ (En np n V n f 

" UW EnPnV n p n V n f ' 



(17) 



The uncertainty in our knowledge of the number of pho- 
tons is the square root of this variance. Averaging this 
uncertainty over all trajectories therefore tells us, on av- 
erage, how accurately we will have determined the num- 
ber of photons at a later time. To calculate the value 
of the uncertainty for each trajectory, averaged over all 
trajectories we must multiply a n (t) w by the probability 
for each final state and average over all the final states. 
The probability measure for the final states, p{t), is given 
by the Wiener measure multiplied by the norm of the fi- 
nal state, Tr{V(t)p(0)Vt(i)}. This probability measure 
is not in general Gaussian in W(t), but a weighted sum 
of Gaussians, one for each n. Performing the multipli- 
cation, we obtain the average conditional uncertainty in 
photon number as 




in which 



V„ 

dP w 



n)p n pmV n V m dP w , (18) 



(19) 
(20) 



-4ktn' 1 +2V2knW 



-W 2 /(2t) dw _ 



We note that (a n (t) w ) may be written as a function of 
t = kt, being the time scaled by the measurement con- 
stant. Hence, as we expect, increasing the measurement 
time has the same effect on (cr n (t) w ) as increasing the 
measurement constant. We evaluate (cr n (r) w ) numeri- 
cally for an initial thermal state, and an initial coherent 
state, and display the results in figure [1] We have chosen 
the initial states so that they have the same uncertainty 
in photon number, with the result that the mean num- 
ber of photons in each of the two states is quite different. 
The results show the decrease in uncertainty with time, 
which is seen to be only weakly dependent upon the ini- 
tial state. 
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Conditional Uncertainty for a QND Measurement of n 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



FIG. 1: The conditional uncertainty in photon number av- 
eraged over all trajectories, {a„(t) w }, is plotted here against 
the dimensionless scaled time, r = kt. The dotted line cor- 
responds to an initial coherent state, and the solid line to an 
initial thermal state. Both initial states were chosen to have 
= 20, giving the thermal state a mean photon number 
of (n) = 4, and the coherent state a mean photon number of 
(n) = 20. The photon number distributions for the two initial 
states are displayed in the inset. 



C. A measurement of momentum in a linear 
potential 

The simplest mathematically non-trivial case occurs 
when the commutator between A and B, while non-zero, 
is such that it commutes with both A and B. A phys- 
ical situation to which this corresponds is a continuous 
measurement of the momentum of a particle in a linear 
potential. If we denote the position and momentum op- 
erators for the particle as Q and P respectively, then the 
Hamiltonian is given by 



H = ^-P 2 - FQ, 
2m 



(21) 



in which m is the mass of the particle and F is the force 
on the particle from the linear potential. In this case we 
have 



A 
B 



2hm 
2fcP, 



2k)P 2 



iF 



(22) 
(23) 



in which k is again the measurement constant. 

Returning to Eq. (fTU|) we see that to obtain a solution 
we must pass all the exponentials containing the operator 
B to the right through the exponentials containing the 
operator A. In order to perform this operation we need 
a relation of the form 



e B e A 



e A e c . 



(24) 



For the present case the required relation is simply given 
by the Baker-Campbell-Hausdorff formula [l], [HJ 



B A A B -\A,B] 

e e = e e e 1 J . 



(25) 



Using this relation to propagate successively all of the 
exponentials containing B to the right in the product in 
Eq. (fT0|) we obtain 



JV 

IK 



e AAt e BAW„ 



) = cxp [AN At] exp 



N 



Bj2&W n 



x cxp 



N 



-[A,B]AtJ2(n-l)AW n 



71=1 



(26) 



All that remains is to calculate the joint probability den- 
sity for the random variables. The first is the Wiener 
process, and the second is 



dW(t'){27) 



Y(t) = lim At V(n - 1)AW„ = / t' 

At ^° n=l J 

Clearly these are both Gaussian distributed with zero 
mean and all that we require is to calculate the covari- 
ances (Y(t) 2 ) and (W(t)Y(t)) to determine completely 
the joint density at time t. Using (AW n AW n ) = 8 nm At 
these quantities are easily obtained: 

N 

(Y(t) 2 ) = lim o Ai]T((n-l)At) 2 At 

^ n=l 

t' 2 dt' = t 3 /'S, (28) 

JV 

lim V((n- I) At) At 



At^o 

rt 



(W(t)Y(t)) 



t' dt' = t 2 /2. 



(29) 



The state at time t, under the evolution described by the 
stochastic equation, is therefore 



\^{t)) w = e At e BW ^e^ A ^ Y ^\m) ) 



(30) 



where the joint probability density for W and Y at time 
t is given by 



P V (W,Y) 



12 



27rf 2 



exp 



2 9 6 9 6 
— W 2 - -^Y 2 + ^rWY 

t t 3 t 2 



Note that to obtain the probability density for the final 
state, this must be multiplied by the norm of the state 
at time t. 

Returning to the specific case of a particle in a linear 
potential, we may now obtain results for various quanti- 
ties of interest. Writing the evolution operator in terms 
of the momentum and position operators we have 



= exp 
x exp 



(( 



2k)P 2 



iF 



Q)t 



2hm h 
2k(PW(t) + FY(t)) |-0(O)>. (31) 
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Using the Zassenhaus formula [23j to disentangle the ar- 
gument of the first exponential we may rewrite this in 
the more convenient form 



exp 



x exp 



iF 



-Qt 



exp [ri{-P 2 t - PFt 2 - F 2 i 3 /3)] 



2k(PW(t) + FY(t)) |^(0)), 



(32) 



in which ij = (at + ^or those not familiar with 
the Zassenhaus formula, it is complementary to the BCH 
formula. While the BCH formula shows how to write an 
exponential of the sum of two operators as a product 
of exponentials of the operators and their commutator 
(or in more complicated cases repeated commutators of 
the two operators), the Zassenhaus formula shows how 
to write the product of exponentials of two operators as 
the exponential of the sum of the operators and repeated 
commutators. 

Let us first take an arbitrary initial state, writing it in 
the momentum eigenbasis so that we have 



Mo)) 



*(p)b) dp 



\*(p)\ 2 dp=l. (33) 



Using Eqs. (|32|) and (|33|) in Eq.© to calculate the mo- 
ments of p given by first averaging over all trajectories 
(that is, the moments which would be given by the equiv- 
alent master equation) we readily obtain 



( P (tT) = ((p(o) + Fty 



(34) 



In particular, for any initial state, \ip(0)}, the average 
value of the momentum at time t, (p(t)}, is simply shifted 
from the initial value by the impulse Ft. The variance 
of the momentum at time t, a 2 (t) = (p(t) 2 ) — (p(t)) 2 , 
remains equal to its original value. That is, the uncer- 
tainty introduced into the position of the particle by the 
momentum measurement does not feed back into the mo- 
mentum, even though the momentum does not commute 
with the Hamiltonian. This is because while the momen- 
tum determines the position at a later time, the converse 
is not true. These results for the moments are easily 
chocked using the equivalent master equation. 

Now let us consider the conditional variance of the 
momentum at time t averaged over all trajectories. In 
the previous section we calculated the conditional un- 
certainty, being the square root of the variance, and av- 
eraged this over all trajectories. Here however, we will 
find that the conditional variance is independent of the 
trajectory taken, and depends only on the measurement 
time. This will also be true of the example which we 
will treat in the next section. In this case clearly it does 
not matter if we first average the conditional variance 
over the trajectories, and then take the square root, or if 
instead we average the conditional uncertainty, because 
the averaging procedure is redundant. However, in gen- 
eral the two procedures are not equivalent. We will de- 
note the conditional variance by (a 2 (t) w ). As the uncer- 
tainty in position does not feed back into the momentum, 



we expect that this variance should steadily decrease to 
zero. This is because during a trajectory our knowledge 
of the momentum steadily increases so that the distribu- 
tion over momentum becomes increasingly narrow. To 
perform this calculation we take the initial state to be the 
minimum uncertainty wave-packet given by the ground 
state of a harmonic oscillator of frequency to. The aver- 
age values of the position and momentum of the particle 
are both zero in this state and the respective variances 
are 



(Q 2 ) 



. ( p ) - 



*p(0), 



and in momentum space the state may be written 



|^(0)) 



1 



ihuj 



e 

— oo 



-P 2 /(2mM|p) dp. (35) 



The moments of momentum for each trajectory are given 
by 



m)\p n m))^ 
m)W)u 



(36) 



and we calculate the first and second to give <?p(t) y 
(p(i) 2 )„ - (p{t)) 2 w . We obtain 



l fl Jw] " l + 8fco2(0)f 



(37) 



This is independent of W and Y and hence independent 
of the trajectory. It is therefore unnecessary to average 
over the final states. Indeed (cr 2 (i) w ) decreases steadily 
from the initial value to zero as t — > oo as we expect from 
the discussion above. This means that while the average 
value of momentum is determined by the measurement 
results, the error in our estimate of the momentum at 
time t is not. 



D. A quadrature measurement with a general 
quadratic Hamiltonian 

We now consider an LSE in which the commutator 
[A, B] does not commute with either A or B. As in 
the previous example, let P and Q be, respectively, the 
canonical momentum and position operators for a single 
particle so that they obey the canonical commutation re- 
lation [Q,P] = iH. With this definition we will take A 
and B to have the following forms: 

A = aP 2 + jQ 2 + £QP + r/P + £Q, (38) 
B = kQ + nP, (39) 

where a, 7,77, £, k and k are complex numbers. This ex- 
ample applies to an optical mode of the electromagnetic 
field, including classical driving and/or classically driven 
subharmonic generation [2~i| and for which an arbitrary 
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quadrature is continuously measured. It also applies to 
the situation of a single particle, which may feel a lin- 
ear and/or harmonic potential, and which is subjected 
to continuous observation of an arbitrary linear combi- 
nation of its position (25| and momentum. 

To obtain an evolution operator for the LSE with this 
choice of the operators A and B, we require, as before, 
a relation of the form given by Eg . (|2"4")1 . To derive this 
relation we proceed in the following manner. 

First we may use the Baker-Campbell-Hausdorff ex- 
pansion [![, or alternatively solve the equations of mo- 
tion given by dB/de = [A,B], to obtain an expression 
for e eA Be~ eA . The result is 



e- eA eBe eA 



ef 1 (e)Q + ef 2 (e)P + ef 3 (e), (40) 



in which 



h(e) = j(-2K>y + k£)S + kC, 
f 2 (e) = i(2ka-K0S + KC, 
h(e) = i(A^ + 2ftaC-«Ce-2«7Tj)[C-l] 



A 2 

+~(kr) + K()S. 
A 



(41) 
(42) 

(43) 



In these expressions C = cosh(iSAe), S = sinh(ifiAe) and 

f(e- £A eBe eA ), (44) 



A = \Ji 2 — 4a-f. Using the relation 



e-* A f(eB)e 
we obtain from Ea.(|40p 



-tA n eB n eA 



ef 1 (e)Q+ef 2 (e)P f ,ef 3 (e) 



(45) 



Multiplying both sides of this equation on the left by e 



<A 



eB „nA 



f-A a eD(t) 



e e 



we obtain a relation of the form e e 
we require. 

We see from the above procedure that the relation in 
Eq. (f2"4"| may be obtained so long as a closed form can be 
found for the solution to the operator differential equa- 
tion dB/de = [A,B]. Clearly this is straightforward if 
this equation is linear in B, which is true in the example 
we have treated here, and is sometimes possible in cases 
in which the equations are non-linear. 

In addition, for this example we also require the BCH 
relation in the form 



£ A e B = e A+B+\[A,B}^ 



(46) 



This is so that we can sum up in one exponential the 
operators that result from swapping e Aw ^ B and e nAtA . 

Using the expressions derived above, with the replace- 
ments e = n At and e = A W„ , for each n from 1 to N, by 
repeatedly swapping the exponentials containing B with 
those containing A as in the previous example, we obtain 



N 



lim TT(< 



AAt BAW n 



e At e X 1 (t)Q+X 2 (t)P 



xe x 3 (t) e ihZ(t) 



in which the classical stochastic variables Xi and Z, are 
given by 



(48) 



Z(t) = I ' h{t')X 2 {t')dW{t')- [ hWXtWdWtf), 
Jo Jo 

where the expressions for the fi are given above, and 
the integrals are Ito integrals. The Xi arc Gaussian dis- 
tributed with zero mean, and their covariances are easily 
calculated as in the previous example: 

{X i {t)X j {t)) = (f i {t')f j {t')dt l . (49) 

In addition, the two-time correlation functions for these 
variables are also easily obtained analytically. In partic- 
ular we have 



(X^X^t)) = 



min(t,r) 



(47) 



/i(f)/i(*') dt'. (50) 



However, Z(t) is not Gaussian distributed. We are not 
aware of an analytic expression for this variable, so that 
its probability density may have to be obtained numer- 
ically. We note in passing, however, that in some cases 
double stochastic integrals of this kind may be written 
explicitly in terms of products of Gaussian variables 0] • 
We note also that Z determines only the normalisation of 
the final state, and not the state itself. The normalised 
state at time t is therefore independent of Z, and we 
examine the consequences of this in appendix [5] 
We may now write the state at time t as 

|^(i)) w = eM^m+x^p^+mzM^y (51) 

Hence even though values for averages over all trajec- 
tories may in general have to be calculated numerically, 
the evolution operator provides us with information re- 
garding the type of states that will occur at time t. In 
particular, if the initial state is Gaussian in position (and 
therefore also Gaussian in momentum), then as each of 
the exponential operators in the above equation trans- 
form Gaussian states to Gaussian states, we see that the 
state of the system remains Gaussian at all times. The 
mean of the Gaussian in both position and momentum 
change with time in a random way determined by the 
values of the stochastic variables. 

We will shortly consider a particular example: that of a 
harmonic oscillator undergoing a continuous observation 
of position, and use this evolution operator to calculate 
the conditional variance for the position at time t. We 
will take the initial state to be a coherent state, which 
is a Gaussian wave packet. This conditional variance 
does not depend upon the trajectory, but simply upon 
the initial state and the measurement time, as indeed we 
found to be the case for the momentum measurement in 
section III CI 
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Let us first show that for an initial coherent state the 
conditional variance of any linear combination of position 
and momentum is independent of the trajectory for all of 
the cases covered by the evolution operator in Eq. (f5"Tj) . 
To do this we must calculate the effect of this evolution 
operator on a coherent state. Clearly the effect of the 
right-most exponential operator is at most to change the 
normalisation, which effects neither the average values 
of position and momentum, nor the respective variances. 
The effect of the next exponential, being linear in P and 
Q, is calculated in appendix [Bl We find that it changes 
the mean values of the position and momentum, and al- 
ters the normalisation, but the state remains coherent 
in that the position variance (and hence the momentum 
variance) is unchanged. Finally, the effect of the expo- 
nential quadratic in P and Q is calculated in appendix IB1 
We find that this operator modifies the position variance. 
However, as the operator does not contain any stochastic 
variables, and as the manner in which it changes the po- 
sition variance is independent of the mean position and 
momentum, we obtain the result that the effect on the 
position variance, and hence the variance of any linear 
combination of position and momentum, is trajectory in- 
dependent. 

Let us now consider a harmonic oscillator in which 
the position is continuously observed. This situation 
has been analysed by Belavkin and Staszewski using the 
equivalent non-linear equations [l3|. The operators A 
and B in this case are given by 



*-(—) 

\2Hm J 



D 



p2 



2h 



2kQ, 



2k ) Q\ (52) 
(53) 



in which m is the mass of the particle, u is the frequency 
of the harmonic oscillation, and k is the measurement 
constant for the continuous observation of position. Tak- 
ing the initial state to be coherent, and denoting it \a), 
the initial position wave-function is given by 



After some algebra this may be written as 

iztanh(zujt) — 1 
tanh(zLut) — iz 



n 2 
s = s iz- 



(58) 



in agreement with that derived by Belavkin and 
Staszewski. The conditional variance for x at time t is 
given by 



1 



4Re[s' 2 ] ' 



(59) 



As t tends to infinity, Eq. (j58]) gives a steady state value 
for the conditional variance, which is 



1 



4Im[z] 



1 + 1 



(60) 



The parameter r is a dimcnsionlcss quantity which gives 
essentially the ratio between the frequency of the har- 
monic oscillator, and the rate of the position measure- 
ment. We may view the dynamics of the position vari- 
ance as being the result of two competing effects. One 
is the action of the measurement which is continuously 
narrowing the distribution in position, and consequently 
widening the distribution in momentum. The other is 
the action of the harmonic motion, which rotates the 
state in phase space, so converting the widened momen- 
tum distribution into position. Depending on the relative 
strengths of these two processes, determined by the di- 
mensionless constant r, a steady state is reached in which 
they balance. If the rate of the measurement is very fast 
compared to the frequency of the oscillation (correspond- 
ing to r <C 1), then the localisation in position is much 
greater than it would be for an unmonitored oscillator, 
and in that case we succeed effectively in tracking the 
position of the particle. However, if the frequency of os- 
cillation is much greater than the rate of localisation do to 
the measurement, then the steady state position variance 
remains essentially that of the unmonitored oscillator. 



ma) 



2s2 

7T 



1/4 



-s 2 x 2 +2sxa-±(\a\ 2 +a 2 ) 



(54) 



where s 2 = mw/(2fi). Using the results in appendix IBl 
we find that the coefficient of x 2 at a later time t is given 
by 



1 - 21 



3-21 



1 



1-21 
1 + 21 



where 



-1/2 



rz coth(zcut) + (1 + ir) ' 
and we have defined the parameters 

™, ,2 



1 



2hk 



(55) 



(56) 



(57) 



III. CONCLUSION 

We have presented a method for obtaining evolution 
operators for various classes of stochastic equations de- 
scribing linear quantum trajectories, and applied this to 
a number of physical examples pertaining to physical 
systems subjected to the continuous projection measure- 
ment of an observable. We have shown how the com- 
plexity of the stochastic equations governing the random 
variables which appear in the evolution operator depends 
upon the commutation relations between the operators 
appearing in the LSE. For the case in which both these 
operators commute with their commutator, probability 
densities for the random variables may be obtained an- 
alytically. We have also shown that in cases in which 
the commutation relations are more complex it is some- 
times still possible to obtain an explicit evolution opera- 
tor. This is possible even in cases in which the classical 
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stochastic integrals, or cquivalently the stochastic equa- 
tions, governing the random variables which appear in 
this operator are too complex to solve analytically. 
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APPENDIX A: ELIMINATING VARIABLES 
WHICH AFFECT ONLY THE FINAL 
NORMALISATION 



We found in section III Dl that not all the random vari- 
ables which appear in the evolution operator are Gaus- 
sian distributed. This result is surprising because it has 
been shown previously, using the non- linear equations, 
that for an initial Gaussian state, the probability density 
for the conditional mean position and momentum, and 
therefore for the final state, are Gaussian distributed for 
this case [l2| . These two results may be reconciled due to 
the fact that the non-Gaussian variable in the evolution 
operator given in Eq. (|5ip affects purely the normalisation 
of the final state, rather than the state itself. 

Let us assume that we have an initial state \ip), and 
an evolution operator which is a function of the random 
variables X and Z (which may in general be vector val- 
ued). We let the random variable Z determine only the 
normalisation of the final state, so that the evolution op- 
erator may be written 



V(X,Z,t)=0(X,t)f(Z,t), 



(Al) 



where O is an operator valued function, and / is simply 
a complex valued function. The unnormalised state at 
time t is then given by 



m))™ = o(x,t)f(z,tM. 



(A2) 



Clearly once we have normalised that state at time t, it 
is no-longer dependent upon Z. In particular the nor- 
malised state is given by 



\m) w = 



o(x,tm, 



(A3) 



y/^\O^X,t)0(X,tM w 

The probability density for the final state is 

P(X,Z,t) = (i>(t)\i>(t)} w P w (X,Z,t), (A4) 

in which P W (X, Z,t) is the probability density given by 
the Wiener measure for the variables X and Z . However, 



seeing as the normalised state depends only upon X, we 
require for all calculations only the marginal probability 
density for X. Denoting this marginal density also by P, 
we have 



P{X,t) = J P(X, Z,t) dZ. 



(A5) 



In certain cases the probability measure for the nor- 
malised state may therefore be Gaussian, even though 
the measure for the unnormalised state is not. However, 
as P(X,Z,t) contains a factor of the norm of |-0(i)) w , 
the probability measure for the output process will, in 
general, only be Gaussian if the norm is Gaussian in X. 
Clearly the norm is Gaussian in X for initial Gaussian 
states in the case we investigate in section ITlDl 



APPENDIX B: THE EFFECT ON A COHERENT 
STATE OF EXPONENTIALS LINEAR AND 
QUADRATIC IN P AND Q 

We first calculate the effect of an operator of the form 
e »P+K.Q ( B1 ) 

on a coherent state |a). The coherent state is defined as 
the eigenstate of the annihilation operator a, such that 



a\a) = a\a), 



and 



Imuj I 1 
a = \ I -ri- x + i\l — p. 



2h 



2hmuj 



(B2) 



(B3) 



Here m and tu are the mass and frequency of a harmonic 
oscillator which serves for the purposes of defining the 
coherent state. In particular we are interested in the 
position wave- function of the result. We therefore wish 
to calculate 



<a#> = (x\e vP+,lQ \a), 



(B4) 



where \x) is an eigenstate of the position operator Q such 
that 



Q\x) = x\x). 



(B5) 



Note that in general \4>) will not be normalised. To per- 
form this calculation we will need the BCH formula given 
in Eq. (|46p . and the position wavefunction for a coherent 
state, 



(x\a) = 

7T 



2 S ^ 1/4 



-s 2 x 2 +2sxa-±(\a\ 2 +a 2 ^ 



-s x -\-2sxa — a —ia- r oti 



(B6) 
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where 



in which 



s 
a 



imu 
~2h 

a r + ioli 



(B7) 
(B8) 



Note that this expression contains the phase factor 
— ia r on. This is left out in many texts, but is essen- 
tial for consistency with the completeness relations for 
the position states. We also require the inner product of 
two coherent states, 



(B9) 



and the well known integral formula 

, Re [a] > 0. (BIO) 



e - ax --p x dx = d /^ p /3 2 /(4a) 



We proceed first by rewriting the exponential in terms of 
annihilation and creation operators, so that we have 



a uP+fiQ _ fla+0a+ 



e e 



in which 



2muj 



mhuj 



(Bll) 



(B12) 



We may now use the completeness relation for the coher- 
ent states to obtain 



(x|V>) = (x\e 



= (x\a + one 



|e^V»e^ 2 d 2 /? 



9q+0/3*+90/2 ^2^ 



l \<t>\' 2 +Re[a<t>*]+0a+9ct>/2 



(B13) 



We see that the state remains coherent, although it is no 
longer normalised, and is shifted in phase space by <j). 

We now wish to calculate the effect of an operator of 
the form 



D vP 2 +CQ 2 +iQP 



(B14) 



on a coherent state. This time we require to calculate 

(a#) = (x|e" p2+CQ2+ « QP |a). (B15) 

For this calculation we will need the disentangling the- 
orem for the exponential of a general quadratic form of 
the annihilation and creation operators, which is given 
by M 



ua 2 -\-va^ 2 -\-wa^ a (x^^)/2 Za^ 2 %a^a ma 2 



(B16) 



/coth(/) - w 



In 



/ 



/cosh(/) - wsinh(/) 
u 



m 



/coth(/) - w 
\/ w 2 — 4uv. 



(B17) 

(B18) 
(B19) 
(B20) 



First of all rewriting the exponential containing P and Q 
as an exponential in the annihilation and creation oper- 
ators, we have 



(x\e r - 



in which 



» = 



( Qi -qmhw £h 
= v = i — 



V Imuj 



r/mhui I 



\a) (B21) 

(B22) 
(B23) 



We now proceed by using the disentangling theorem, 
and employing the completeness relation for the coherent 
states. 

(x\1i) = {x\e {w+x)/2 e la ' 2 e xa ^ a e ma2 \a) 

= \ JJ(x\l3}(f3\e^ w+ ^/ 2 e la ' 2 e xa " a e ma2 \a) d 2 /3 

= ~ J J (x\/3)((]\ae x ) e 10 ' 2 d 2 f3 

x e i\a\ 2 (\e 2x \-i)+m a 2 e ( x -w)/2 _ (B24) 

Performing the integral over the real and imaginary parts 
of a, we obtain 



1 



2 \ 4 



2s 



vT+27 V ^ 



x exp < —s 2 x 2 



x exp < 2sxae x 



e 2 -|a| -ma & (x-w)/2 



1 - 21 



x exp < a 2 e 2x 



3-2/ 
1 

3~~2l 
1 



1 + 2- 



1 - 21 



1 + 2- 



1 + 21 
1 - 21 



3-21 



1 + 21 
2 



1 



2 1 + 2/ 



(B25) 



It is easily verified that this reduces to (x\a) as required 
when we set / = x = m = 0. 
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